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Abstract. The differential event rate for direct detection of dark matter, both the time 
' (— I ' averaged and the modulated one due to the motion of the Earth, are discussed. The 

Q | calculations focus on relatively light cold dark matter candidates (WIMP) and low energy 

I transfers. It is shown that for some WIMP masses the modulation amplitude may change 

sign. This effect can be exploited to yield information about the mass of the dark matter 
~h candidate. 

>: 

$5 " 1- Introduction 

The combined MAXIMA-1 p], BOOMERANG [2], DASI [3] and COBE/DMR Cosmic 
lO \ Microwave Background (CMB) observations [4J imply that the Universe is flat [5] and that 

most of the matter in the Universe is Dark [6j[7], i.e. exotic. Combining the data of these 
. quite precise experiments one finds: 

(N 

n b = 0.0456 ± 0.0015, n cnM = 0.228 ±0.013, Q A = 0.726 ± 0.015. 
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Since any "invisible" non exotic component cannot possibly exceed 40% of the above 
fi CDM [8], exotic (non baryonic) matter is required and there is room for cold dark matter 
candidates or WIMPs (Weakly Interacting Massive Particles). 

Even though there exists firm indirect evidence for a halo of dark matter in galaxies 
from the observed rotational curves, see e.g the review [9], it is essential to directly detect 
such matter. The possibility of such detection, however, depends on the nature of the dark 
matter constituents and their interactions. 

Since the WIMP's are expected to be extremely non relativistic, with average kinetic 
energy (T) « 50 keV(m WIMP /100 GeV), they are not likely to excite the nucleus. So they 
can be directly detected mainly via the recoiling of a nucleus (A,Z) in elastic scattering. 
The event rate for such a process can be computed from the following ingredients [10] : 
i) The elementary nucleon cross section, ii) knowledge of the relevant nuclear matrix 
elements 1 12 j . obtained with as reliable as possible many body nuclear wave functions, 
hi) knowledge of the WIMP density in our vicinity and its velocity distribution. 

The nucleon cross sections can also be extracted from the data of event rates, if and 
when such data become available. From limits on the event rates, one can obtain exclusion 
plots on the nucleon cross sections as functions of the WIMP mass. In the standard nuclear 
recoil experiments, first proposed more than 30 years ago [13] , one has to face the problem 
that the reaction of interest does not have a characteristic feature to distinguish it from 
the background. So for the expected low counting rates the background is a formidable 



problem. Some special features of the WIMP-nuclear interaction can be exploited to reduce 
the background problems. Such are: 

i) the modulation effect: this yields a periodic signal due to the motion of the earth 
around the sun. Unfortunately this effect, also proposed a long time ago P3] and 
subsequently studied by many authors HQ [ISl H3 HHl QSl I2D1 I2H 1221 123], is small and 
becomes even smaller than 2% due to cancelations arising from nuclear physics effects, 

ii) backward-forward asymmetry expected in directional experiments, i.e. experiments 
in which the direction of the recoiling nucleus is also observed. Such an asymmetry has also 
been predicted a long time ago [M] , but it has not been exploited, since such experiments 
have been considered very difficult to perform, but they now appear to be feasible[10j. iii) 
transitions to excited states: in this case one need not measure nuclear recoils, but the 
de-excitation 7 rays. This can happen only in very special cases since the average WIMP 
energy is too low to excite the nucleus. It has, however, been found that in the special 
case of the target 127 I such a process is feasible [25] with branching ratios around 5%, (iv) 
detection of electrons produced during the WIMP-nucleus collision [261 E7J an d v) detection 
of hard X-rays produced when the inner shell holes are filled [28J. 

In the present paper we will limit our attention to the standaerd recoil experiments and 
study the differential event rates, both time averaged and modulated, in the region of low 
energy transfers, as in the DAMA experiment |29[ 130] , focusing our attention on relatively 
light WIMPS [311 E21 E3]. Such light WIMPs can be accommodated in some SUSY models 
[34] . We will also present some results on the total rates as well. We will employ here the 
standard Maxwell-Boltzmann (M-B) distribution for the WIMPs of our galaxy and we will 
not be concerned with other distributions [35} 1361 [371 EH] , even though some of them may 
affect the modulation. The latter will be studied elsewhere. We will explicitly show that the 
modulation amplitude, entering both the differential and the total rates, changes sign for 
certain reduced WIMP-nuclear masses. As a result such an effect, if and when the needed 
data become available, may be exploited to infer the WIMP mass. 



2. The formalism for the WIMP-nucleus differential event rate 

This formalism adopted in this work is well known (see e.g. the recent reviews [391 140]). So 
we will briefly discuss its essential elements here. The differential event rate can be cast in 
the form: 

dR. dR , dH. 

dQ U = -dQ U+ dQ lACOSa (1) 

where the first term represents the time averaged (non modulated) differential event rate, 
while the second gives the time dependent (modulated) one due to the motion of the Earth 
(see below). Furthermore 

— \ A = Px mt J ( ^ V V< v 9 ' > A 2 1 dt 
dQ m x Am p n \^ p J Qo(A) du 

™ U = ft^fit)"^^^ 1 * (2) 
dQ m x Am p \^pj Qo(A) du 

with with fjL r (/i p ) the WIMP-nucleus (nucleon) reduced mass, A is the nuclear mass number 
and a n is the elementary WIMP-nucleon cross section. m x is the WIMP mass and mt the 
mass of the target. Furthermore one can show that 

|- = y|a 2 F 2 (u)^o(av^), fj = y|a 2 F 2 («)* 1 (aV^ (3) 



with a = (v^r^o) -1 ) Vq the velocity of the sun around the center of the galaxy and b the 
nuclear harmonic oscillator size parameter characterizing the nuclear wave function, u is 
the energy transfer Q in dimensionless units given by 

u = TTTIT ' Qo(A) = [m p Ab 2 ]~ l = 40vl" 4/3 MeV (4) 

and F(u) is the nuclear form factor. Note that the parameter a depends both on the WIMP 
, the target and the velocity distribution. Note also that for a given energy transfer Q the 
quantity u depends on A. 

The functions ^o(ai/u) and ^>i{a^/u) for a M-B distribution take the following form: 

*o0) = ^(erf(l -x) + erf(x + 1) + «fc(l - y csc ) + erfc(y csc + 1) - 2) (5) 

^i(x) = -6 (-erf(l - x) - erf(x + 1) - erfc(l - y csc ) - erfc(y csc + 1) 

2 e -( x ~ 1 ? 2e~( x+1 ) 2 2e-^ csc_1 ) 2 2e~( J/csc+1 ) 2 \ 

+ 7^ + 7= 7= 7= + 2 (6) 

vvr xPk x Fk xPk 



where erf(x) and erfc(x) are the error function and its complement respectively, 5 ~ 0.135 
is the ratio of the velocity of the Earth around the sun to that of the sun around the center 
of the galaxy and a is the phase of the Earth (a = 0, around June 3nd). 
Sometimes we will write the differential rate as: 

= ( -Y^^^-L- ( ±(l+H(ay/u)casa) (7) 

dQ 1 m x Am p V/W Qo(A) \du y v ' J w 

In this formulation H(a^/u), the ratio of the modulated to the non modulated differential 
rate, gives the relative differential modulation amplitude. 

The function H{a^/u) is shown in Fig. [TJ It is independent of the nuclear physics and 
depends only on the reduced mass and the velocity distribution. They are thus the same 
for both the coherent and the spin mode. Note that H(a^Ju) can take both positive and 
negative values, which affects the location of the maximum. 



3. Some results on differential rates 

We will apply the above formalism in the case of Nal, a target used in the DAMA experiment 
[29| [30] . The results for the Xe target are similar [31] . The differential rates \ a and 

^tjU) for each component (A = 127 and A = 23) are exhibited in Fig. [2][3l Following 
the practice of the DAMA experiment we express the energy transfer is in keVee using the 
phenomenological quenching factor [41} 142] . The nuclear form factor has been included (for 
a heavy target, like 127 I or 131 Xe, its effect is sizable even for an energy transfer [10] of 10 
keV). The differential rate for the spin mode for low energy transfers is similar to those 
exhibited in Figs [2][3l since the spin form factors are similar. They are, of course, simply 
scaled down by A 2 , if one takes the spin cross section, a combination of the nuclear spin ME 
and the nucleon spin amplitudes, to be the same with the coherent nucleon cross section, 
i.e. <7^™ ear = 10 _7 pb. For the actual spin nucleon cross sections extracted from experiment 
see [S] and [H S3 M\ ■ 

The functions H(ay/u) cos a for each target component are shown in Figs [U- [5] as a 
function of a for various low energy transfers. The corresponding quantities for the spin 




Q ^keV 

Figure 1: The function H(ay/u) entering the modulated differential rate as a function of 
the recoil energy for a heavy target, e.g. 127 I. Note that this is independent of the form 
factor. The solid, dotted, dot-dashed, dashed, long dashed and thick solid lines correspond 
to 5, 7, 10, 20, 50 and 100 GeV WIMP masses. 
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Figure 2: The differential rate ^j, as a function of the recoil energy for a heavy target, e.g. 
127 I (a) and the amplitude for the modulated differential rate ^ty (b), assuming a nucleon 
cross section of 10 pb. The solid, dotted, dot-dashed, dashed, long dashed and thick solid 
lines correspond to 5, 7, 10, 20, 50 and 100 GeV WIMP masses. Note that is given in 
absolute units. 
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Figure 4: The modulation H(ay/u) cos a with an energy transfer of 2 keVee (a) and 5 keVee 
(b) for a heavy target (I or Xe). The solid, dotted, dot-dashed, dashed, long dashed and 
thick solid lines correspond to 5, 7, 10, 20, 50 and 100 GeV WIMP masses. Note that 
for some wimp masses on June 2nd the amplitude becomes negative (location of minimum 
rate). Note that the modulation is given relative to the time averaged rate. 



mode are almost identical. We see that for certain values of the WIMP mass the modulation 
amplitude changes sign. This may perhaps by exploited to extract information on the WIMP 
mass from the data. A similar behavior has been found by considering various halo models 
and different minimum WIMP velocities [22\ 123] . Sometimes, as is the case for the DAMA 
experiment, the target has many components. In such cases the above formalism can be 
applied as follows: 

— — \a — > y^-^i-r^U-; u ^ Ui, X% = the fraction of the component Aj in the target 
dQ ' dQ 

i 

(8) 

The function H(a,y/u) cos a for Nal as a function of a is exhibited in Fig. [6j 
4. Some results on total rates 

For completeness and comparison we will briefly present our results on the total rates. 
Integrating the differential rates discussed in the previous section we obtain the total time 
averaged rate Rq, the total modulated rate H and the relative modulation amplitude/i given 




Figure 5: The same as in Fig. [J] for a light target (Na or F). 




Table 1: Some total event rates for some special WIMP masses and energy thresholds. The 
coherent nucleon cross section of a n = 10 pb was employed. 





m WIMP 


«o(I) 


HP) 


hp) 




fi(Na) 


h(Na) 


flo(Nal) 


ff(Nal) 


h(Nal) 


(kcVcc) 


GcV 


kg-y 


kg-y 




kg-y 


kg-y 




kg-y 


kg-y 







SO 


16.3 


-0.311 


-0.019 


1.518 


0.028 


0.019 


14.0 


-0.259 


-0.018 





20 


25.8 


0.285 


0.011 


2.35 


0.050 


0.021 


22.2 


0.249 


0.019 





10 


18.4 


0.356 


0.019 


2.045 


0.046 


0.022 


15.9 


0.309 


0.019 


5 


80 


7.00 


-0.042 


-0.006 


1.133 


0.038 


0.034 


6.11 


-0.030 


-0.005 


5 


20 


2.72 


0.247 


0.091 


1.07 


0.065 


0.060 


2.47 


0.219 


0.089 


5 


10 


0.008 


0.001 


0.187 


0.303 


0.031 


0.103 


0.053 


0.006 


0.114 



by: 

R = Ro + H cos a or R = Rq (1 + h cos a) (9) 

Some special results in the case of low WIMP mass are exhibited in Tables [l][2l From table 
[2] it becomes clear that, for low mass WIMPs, large nucleon cross sections can accommodate 
the data. A similar interpretation holds for the [32] data. In the case of non zero threshold 
one notices the strong dependence of the rime averaged rate on the WIMP mass. Also 
in this case the relative modulation h substantially increases, the difference between the 
maximum and the minimum can reach 20%. This however occurs at the expense of the 
number of counts, since both the time averaged and the time dependent part decrease, but 
the time averaged part decreases faster. So their ratio increases. This can be understood by 
noticing that the cancellation of the negative and positive parts in the differential modulated 
amplitide, becomes less effective in this case. 



Table 2: The same as in table Q] for <r n = 2 x 10~ 4 pb relevant for the DAMA region. One 
sees that, for very low mass WIMPs, large nucleon cross sections are required to obtain the 
rates claimed by the DAMA experiment |30j. 





m WIMP 


Ho (I) 




H(l) 


h(I) 




H(Na) 


h(Na) 


R (NaI) 


if (Nal) 


h(Nal) 


(koVco) 


GcV 


kg-y 




kg-y 




kg-y 


kg-y 




kg-y 


kg-y 







80 


4.07 x 


10 4 


-776 


-0.019 


3.80 X 10 3 


70.2 


0.019 


3.50 X 10* 


-647 


-0.018 





20 


6.43 X 


10 4 


712 


0.011 


5.87 X 10 3 


126 


0.021 


5.54 x 10 4 


622 


0.011 





10 


4.61 X 


10 4 


891 


0.019 


5.11 X 10 3 


115 


0.022 


3.98 X 10 4 


772 


0.019 


5 


80 


1.75 X 


10 4 


-105 


-0.006 


4.83 X 10 3 


95.0 


0.034 


1.53 X 10 4 


-74.6 


-0.005 


5 


20 


6.80 X 


10 3 


617 


0.091 


2.69 X 10 3 


162 


0.060 


6.17 x 10 3 


547 


0.089 


5 


10 


19.4 




3.62 


0.187 


757 


78.1 


0.103 


132 


15.0 


0.114 



5. Discussion 

In the present paper we obtained results on the differential event rates, both modulated and 
time averaged, focusing our attention on small energy transfers and relatively light WIMPS. 
We found that: 

• The relative modulation amplitude crucially depends on the WIMP mass. For small 
masses it exhibits normal behavior, but for large masses it changes sign (minimum in 
June). This effect is more pronounced in the case of heavy targets. We thus suggest 
an analysis of the experiments, in particular the DAMA experiment, along these lines 
to establish the location of the maximum on the a-axis. 

• The relative modulation amplitude depends somewhat on the energy transfer, especially 
at low transfers. 

• For WIMP masses less than 10 GeV, the difference between the maximum and the 
minimum could reach between 20% and 40% for a heavy target, but it is a bit less for 
a light target, depending on the enegy transfer. 

• The relative modulation amplitude for Nal is the weighted average of its two 
components, and in the low energy regime, between 1 and 6 keVee, it does not change 
much with the energy transfer. 

• Once it is established that one actually observes the modulation effect, the sign of the 
modulation may be exploited to infer the WIMP mass. 

For low WIMP mass the total rates depend strongly on the threshold energy, especially for a 
heavy target. The relative modulation in the presence a threshold gets quite large (h ~ 0.2), 
but, unfortunately, this occurs at the expense of the number of counts. It is important to 
compare the relative total modulation in a least one light and one heavy target. For very 
low energy thresholds, if the signs are opposite, one may infer that the WIMP is heavy, 
"t-wimp > 100 GeV. 
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